typeset using JPSJ.sty <ver.l.Ob> 



Thermodynamics and Optical Conductivity of a Dissipative Carrier 

in a Tight Binding Model 

Takeo Kato and Masatoshi Imada 

Institute for Solid State Physics, University of Tokyo, 
7-22-1 Roppongi, Minato-ku, Tokyo 106 

(Received November 5,1997) 

Thermodynamics and transport properties of a dissipative particle in a tight-binding model 
are studied through specific heat and optical conductivity. A weak coupling theory is constituted 
to study the crossover behavior between the low-temperature region and the high-temperature 
region analytically. We found that coherent part around zero frequency in the optical conductiv- 
ity disappears for < s < 2, where s is an exponent of a spectral function of the environment. 
Detailed calculation is performed for ohmic damping (s = 1). In this case, the specific heat 
shows an unusual T-linear behavior at low temperatures, which indicates that the environment 
strongly influences the particle motion, and changes the low-energy states of the dissipative par- 
ticle. The optical conductivity o-(ui) takes a non-Drude form even at zero temperature, and the 
high-frequency side behaves as w ~ 2 , where K is a dimensionless damping strength. The high 
frequency side of the optical conductivity is independent of temperatures, while the low fre- 
quency side depends on the temperature, and behaves as T 2K ~ 2 at high temperatures. We also 
comment on the application of this model to the description of incoherent motion in correlated 
electron systems. 

KEYWORDS: dissipation, specific heat, optical conductivity, ohmic damping, Kubo formula, coherence, incoher- 
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§1. Introduction 

The dissipative dynamics of a quantum particle coupled to a heat bath has been studied for several decades 
in many areas of physics HP The dissipative dynamics has been studied theoretically on a phenomenological 
model in which a particle is coupled to harmonic oscillators. This phenomenological model is called the 
Caldeira-Leggett model. The model has been used to study tunneling effects of dissipative phase motion in 
Josephson junctions.® In the Caldeira-Leggett model, a dissipative particle in a potential V(q) is described 
by the Hamiltonian 
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(1.1) 



(1.2) 



where {qj} and {pj} are sets of positions and momenta of harmonic oscillators. The influence of the envi- 
ronment is determined through the spectral density 

By using .the real-time path integral method, the classical equation of motion of a dissipative particle is 
derived asM 

Mq(t) + M f dt'j(t - t')q(t') + ^-= £(t). (1.4) 



Here, £(t) is a random force, and "/(t) is a damping kernel determined by 



l(t) = Q(*)tt- I"" dio^lcos{ut), (1.5) 
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where 0(i) is a step function. The frequency-independent damping j(t) — jS(t), called the Ohmic damping, 
is obtained by taking the spectral function as J(uj) = M^fuj. 

In the framework of the Caldeira-Leggett theory, various situations have been studied by changing the 
form of the potential. For example, the simplest situation is a case of a flat potential (V(q) = 0). In this 
situation, the particle shows the Brownian motionD'EF The diffusiomxif. a particle has been studied by using 
the real-time path integrals based on the Feynman- Vernon method Jj'Er while the transport properties of a 
dissipative particle can be described by the simple equation of motion 

M(q(t)) + M ( dt'j(t - t')(q(t')) = F(t), (1.6) 



for the average position (q(t)). The random force does not affect the response of the external force F(t), 
because (£(t)) = 0. Particularly for the ohmic damping case j(t) — 7<5(i), the Caldeira-Leggett theory 
corresponds to the naive Drude theory, where 7 -1 is the average collision time. 

On the other hand, the Caldeira-Leggett model has been used to study how the quantum coherence is 
suppressed by dissipation in the double well potential^ which can be truncated to two-state systems called 
the spin-boson model. It has been proved that the spin-boson model shows various behaviors including 
coherent motion at low temperatures for weak damping, and incoherent motion in the limit of high tem- 
peratures and/or strong damping. ThisJeature has been obtained first by an analytical approximation,!!! 
and recently bv numerical simulations aW The spin-boson model with ohmic damping can be related to the 
Kondo modelEf bv the bosonization technique]^ Moreover, the spin-boson_model is relevant to the Fermi 
edge singularities!^ which appears in the X-ray absorption edge anomalytij and the diffusion of a heavy 
particle in metalsO'tSEj'O'H Though these problems have been. solved first by other methods, it has been 
shown that the spin-boson model gives the same result for them.EElEJ Recently, it has been claimed thai, 
spin-boson model with ohmic damping is also relevant to c-axis transport in high-T c superconductors c3 
and a quantum dot coupled to an additional quantum point contact.Er 

In this paper, we consider dissipative dynamics of a quantum particle in the tight-binding model through 
thermodynamics and transport properties. This model for Qhfpic damping is relevant, with certain qualifica- 
tion, to the dissipative motion oLsinall Josephson junctionSjESES' vortices in long Josephson junctionsE^ and 
heavy particles in a metal.E3E3'Ea' We take this model as a phenomenological model to describe the quan- 
tum transport affected by dissipation, which arises from coupling to other degrees of freedop-.j-ln fact, this 
phenomenological introduction is generally justified at least in the high-temperature limit.E3'lHj) We focus 
mainly on the weak-damping region. In this paper, the main aim is to study the crossover behavior from the 
incoherent motion at high temperatures to the quantum motion at low temperatures. At high temperatures, 
the particle motion is described by.the classical master equation in which the tunneling rate is given by the 
same value as two-state systems J^P and the optical conductivity shows an unusual w-dependence because 
of the Fermi edge singularities as seen in two-state systems. On the other hand, in the continuum limit 
which corresponds to the low-temperature and weak-damping limit, we recover the Drude form derived by 
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the equation of motion (1.6) with the mass M replaced by an effective mass M, where M is a curvature 
at the bottom of the particle band. We give the exact formulation for the crossover behavior in the weak- 
coupling region, and connect the two different regions. We also speculate that this study may contribute 
t o un derstand non-Drude forms of the conductivity observed in various strongly correlated systems. (See 

1) 



Here, we introduce the model and review the previous works on this model. In this paper, we consider 
a particle of mass M moving in a periodic potential V(q + a) — V(q) as shown in Fig. 1 (a). Here, a 
barrier height is denoted by Vo, and a frequency of small oscillation around potential minima by u)q. We 
assume Vq 3> huo 3> fc^T, where ks is the Boltzmann constant, and T is a temperature. In this parameter 
region, the system can be reduced to a tight bipd^ng mcjcje|, which has been studied by several authors 



in the framework of the real-time path integral.Ej'EirBBEa' Fisher and Zwerger have studied dynamics 
of a particle in a periodic potential in the weak corrugation limilj-XVo — ► 0), and have related the results 
to the|-tiab|t-binding model through the duality transformation,E3Ea) which is also used in the subsequent 
worksa-EIf In these works, one of the important quantities is the linear mobility fj,i of the particle defined 
in the long-time limit t — > 00 as 

w = l im I im Mil, (1.7) 

F^Ot^oo Ft 
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Fi g. 1 . (a) The periodic potential studied in this paper, and (b) a representative path in the imaginary-time path integral in 



where F is an external force on a particle. The temperature dependence of yn, however, has been exactly 
calculated only for a special value of dimensionless damping strength, K — 1/2, where K = Mja 2 /2-KhE3 l 

In this work, we focus mainly on a weak coupling region K <C 1, and calculate specific heat and optical 
conductivity <j{uj). The linear mobility is also obtained by /i/ = a(ui — ► 0)/e 2 . Thermodynamics of the system 
is formulated by the imaginary-time path integral, and transport properties are studied in the framework of 
the linear response theory. 

The paper is organized as follows. In § |^, we derive exact formal expressions for the partition function 
and the optical conductivity. We also formulate a continuous limit, an incoherent limit and a weak coupling 
theory for arbitrary form of J{oj). For analytical calculations, J(ui) is taken as 

JM^f")" 1 ./^;^), (1.8) 

where 6 S is a dimensionless coupling coefficient, and 5 is a reference frequency. Here, f(uj;uj c ) is a cutoff 
function usually taken as f(uj,uj c ) = e~ w / Wc , while we use a sharp cutoff, f(uj,ui c ) = Q(oj c — uS) in § [Eg] , 
where 0(w) is a step function. In any case, the low-frequency (u> <C w c ) behaviors of the dissipative particle 
are not influenced by the detailed form of the cut-off function. The exponent s in ([l^) mainly determines 
the properties of the heat bath. In § [|, we study the case s = 1 called 'ohmic damping', where the coupling 
constant S s agrees with K. For the ohmic case, application of this calculation to real systems is also discussed 



3.4. In the case s > 1 called 'superohmic damping', and in the case s < 1 called 'subohmic damping', 



calculations are so complicated that we only focus on the continuous limit (§ 2.3) and the incoherent regime 
(§ 2.4). Summary is given in § ||. Throughout this paper, we take h = ks = 1. 



§2. Formulation 



In this section, we give formal expressions for the partition function (§ 2.1), the optical conductivity (§ f^2| 



We also discuss a continuous limit of the tight binding model (§ |2.3|). To deal with this model analytically, 



we focus on two limiting regions: an incoherent tunneling regime (§ 2.4) and a weak damping region (§ [2.5| 



In this section, the derived expressions are applicable to arbitrary form of J(uS). 
2. 1 Partition function 

The exact formal expression for the partition function Z is given in the imaginary-time functional-integral 
represeiilation. The integrals over the heat-bath variables are Gaussian and they can be evaluated ex- 
actly.EJO After integrating out the degrees of freedom for the environment, the partition function is expressed 
as 

Z = Z R fvq( T )exp(-S[q(T)]). (2.1) 

Here, Zr is the partition function of the heat bath, and in following calculation, Zr is removed in order to 
focus only on the damping effects on the particle. The paths q(r) satisfy the periodic boundary condition 
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q(P) = g(0), where (3 — l/T. The effective action S[q(r)] is given by 

S[q(r)] = J dr (±M<f + V(q^j - ± J drdr' 0(r - r')q(r)q(T'). (2.2) 

Damping effects due to the heat bath are described by the last term, which takes a nonlocal form for r. The 
kernel 4>(t) is calculated as 

<t>{r) = - !°° dw^-(D u (0) - D u (t)), (2.3) 
D ^) = \ E ~l^-^ T ( 2 - 4 ) 

71= — OO 

= cosh[q,(/3/2 - |r|)] 

sinh(/3w/2) ' 1 ' j 

where ^„ = 2ixnj (3 is the Matsubara frequency. 

The form of the potential V(q) is shown in Fig. 1 (a). In the parameter region T <C "C Vq, the system 
is effectively mapped to a tight-binding model. In this case, the value of paths q(r) takes q = an for almost 
all r, where n is an integer. A representative path is sJ^WBuan Fig. 1 (b). Though the truncation scheme 
for the paths has already been studied in the literature JiSElliiP we will describe it briefly for self-contained 
description. 

In the dissipationless case, a bare hopping matrix Ao is determined by the instanton action Sq as 

A = Aexp(-S ), (2.6) 

where A is a prefactor of the order of Wq© Hence, we have a tight binding model described by the Hamil- 
tonian 

H = -A Y, (4+iCn + cK+i) . ( 2 -7) 

n 

where 0^(0^) are annihilation (creation) operators of a spinless fermion at site n. The damping effects due 
to the heat bath can be divided into the fast and slow modes as 

J(u) = Ji f (w) + J hf (w) (2.8) 

Ju(w) = J(u)f(u;u c ) (2.9) 

Here, lo c is a cutoff frequency (Aq,T <C lu c "C wo), and f(uj;uj c ) is a cut-off function taken as f(u>;u> c ) <C 1 
for uj c <C w. Then, the fast modes can be treated easily by the Born-Oppenhfiimer approximation. In the 
approximation, the bare matrix element Ao is reduced to a renormalized oneO'E-iP 

*-*>•*(-(; f (210) 

In contrast to the fast modes , th e influence of the slow modes of the heat bath must be studied carefully 
based on the effective action (2.2). In this paper, J\{(lu) is replaced with J(lu), and A and lv c are regarded 
as free parameters (A,T«w c ). 

Because we have already traced out the fast modes which can follow the instanton motion of the time 
scale 1/ojq, each path q(r) is expressed as the sum of the delta functions as 

2m 

q{ T )/a = Y.^ T - Tl )- (2.11) 

1=1 



Here, £j = ±1 specifies the direction of the i-th hopping at r = Tj>. Substituting (2.11) to ( |2.2D , the exact 
formal expression of the partition function is obtained as 



Z = f] A 2 ™ V f dr 2m r m dr 2m _! ■■■ f 2 dn exp 
m=o {(l y Ja Ja Jo 



2m 



.k<l 
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Here, the prime in {£;}' denotes summation in accordance with the constraint 



2»i 



E6 = °- 



(2.12) 



(2.13) 



The expression (2.12) may be interpreted— as the statistical model of the classical particles interacting 
to the r-direction with the potential — c/>(t)& In this point of view, A and £j are regarded as a chemical 
potential and a charge of the classical particle. The inverse temperature f3 = l/T in the original model 
corresponds to the system length in the r direction. At the same time, the constraint (2.13) is interpreted 
as the condition for the charge neutrality. The mapping to the classical statistical mechanics is useful for 
constituting the weak coupling theory. Actually, the 'screening effect' by the charged particles is important 
to describe transport properties of the original model at low temperatures. (See § 2.5.) 



2.2 Optical conductivity _ 
We formulate the optical conductivity o~(u>) in one dimension using the Kubo formulae? 



cr(w) 



2 2 

-e a 



HQ -AM 



i(u + iS) 1 

where 8 is an adiabatic constant, and K, is a kinetic energy 

K = ~ A Y1 ( C n+l C » + C n C «+l) » 
n 

and a current-current correlation function A(w) is defined by 

A(lu) = A(iuj m — » u> + iS). 

r0 



A(iw m ) 



dre l ^(j( T )j(0)), 



Here, the current operator j(r) = c jc r is defined by 



The real part of ct(uj) contains a coherent part expressed by a delta function asol 

Re<r(a;) = DS(u>) + <7 ros (w). 
Here, D is called a Drude weight, and o~ Tes (ui) is a residual incoherent part. From (2.14), we obtain 

D 



, = (-£) - Re A(w->0), 
cr res (a;) = e 2 a 2 



2 ImA(w) 



a; 



(2.14) 
(2.15) 

(2.16) 
(2.17) 

(2.18) 

(2.19) 

(2.20) 
(2.21) 



To evaluate A(w), we consider the imaginary-time correlation function A(r) = 0'( r )j(0)) f° r t > 0. The 
correlation function is formulated by the imaginary-time path integral as 



A(r) = ^ r ^P',(r')ox 1 , [->'</(/.]] 



(2.22) 



Here, the path integral is performed over all the possible paths g(r) with two jumps at r' = 0, t as 

2m 

q(r')/a = ^ 6«J(r' - n) + u5{t') + o'5{t' - r), (2.23) 
i=i 

where 2m + 2 is the number of transitions, and £i,a,a' = ±1 denote the directions of hopping at r' = 77, 0, t, 
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respectively. Then, the path integral is expressed as 

/oo 2m „fi 

W)(--0=E E (-^O^n/ *»(-). (2-24) 

TO=0{f,,<r,o-'}' ' 2 = 1 Jo 

Here, the prime in {£/,<7, a'}' denotes the summation in accordance with the constraint 

2m 

5>+ ff + CT ' = > (2-25) 

z=i 

which comes from the periodic boundary condition for q(r'). Here, we divide A (a;) into two parts as 

A(r)=A 1 (r)-A 2 (r), (2.26) 

where Ai(t ) and ^(t) denote the contribution of paths which satisfy a = — a' and a — a' , respectively. 
From ( |2.22| )-( |2.24| ), we obtain 

2m 



~ . . 2A 2 ^ A 2m ^ fff, 

m =0 {Ci}'- 70 i=l 



' 2m 2m 2m 



exp E telM 7 * _Tfc ) +E " E 6^ T ' - T ) - • (2.27) 



~ 2A 2 ~ A 2m ^ , 



(2m 2m 2m \ 

E - r fe) + E ^ rj ) + E - t ) + ^( t ) 
fc<2 i=l !=1 / 



(2.28) 



Here, the prime in {£/}' and the double prime in {£;}" denote the summations in accordance with ( 2.13| ) 
and the condition 



2 m 



E^+ 2 = 0, (2.29) 



respectively. From (2.27) and (2.28), Ai(r) and A 2 (t) can be interpreted as the partition functions with 
fixed charges at r' = 0, r in terms of the classical partition function (2.12). 

We can calculate a rC s(u) = e 2 a 2 ImA(cj)/cj without the Wick rotation (2.16). We introduce a real-time 
correlation function as 

= -2ImA(T-Mt), (2.30) 
where j(f) = e lHt je~ lHt . Then, the correlation function A(uj) is calculated by 

A(«) - / dtA(t)e iwt . (2.31) 



The above expressions, (2.30) and ( 2.31 ), are used to the calculation in the incoherent region and the 
weak-coupling region. We should note that the causality, A(t) = for t < 0, gives the frequency sum-ruleo 

/oo 
duj Re a{io) = ire 2 a 2 (-£) . (2.32) 
-oo 

The sum rule is always satisfied in the approximations adopted in this paper, since the causality is always 
satisfied. 
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2.3 Continuum limit 

In the dissipationless case, the energy dispersion is given by 



£fc = — 2Acosfca, 



(2.33) 
(A -> oo) 



(2.34) 



where k is a momentum. The lattice model is reduced to the continuum model in the limit a 
by keeping M = l/2Aa 2 constant. In this limit, the energy dispersion is approximated as 

k 2 

e k ~ -2A H 

2M 

Note that the mass M is not necessarily the same as the bare mass M in the original Hamiltonian. In the 
continuum limit A — ► oo, thermal fluctuations may be neglected compared with the band width (T/A — * 0). 
Hence, the continuum limit corresponds to the zero-temperature limit. 

Next, we consider the dissipative case. In the continuum limit, the classical equation of motion for the 
average position (q(t)} is obtained by virtue of Ehrenfcst's theorem as 



(?(*)) + 



dt' 1 (t-t')(q(t')) = Ti E(t), 



(2.35) 



w 
ad 



.e E(t) is an external electric field. The damping kernel 7(4) is determined by the spectral density J(oj) 



lit) 



Ztt 



7(2) = —rj 



duj' 



J(w') 



U!'(U!' 2 + Z 2 ) ' 



(2.36) 
(2.37) 



From the classical equation (2.35), the Fourier transformation of the current j(t) — e(q(t)) gives j(u)) 
a(w)E(u>), and the optical conductivity cr(u>) is given by 

a(u) 



-\UJ-\-8 



M(z + %z)) 

Because of the causality, the sum rule ( 2.32| ) is also satisfied in the continuum limit. 

The above result is valid for any form of the spectral density J(u>). Here, we take J(u>) as 



(2.38) 



J{uj) = Mj s y—J u)Q(lj c - to), 



(2.39) 

where 7 S = 2ir5 s /a 2 Al = 4ir5 s A is a damping frequency, and the cut-off function f(u>;u> c ) in (1.8) is taken 
as the step function 0(w c — u>) for convenience. The characteristic damping frequency 7 is determined 
by 7 = ("/ s LiJ 1 ~' s ) 1 ^ 2 ~ s ' > . To keep 7 constant in the continuum limit A — > 00, the dimensionless coupling 
coefficient 5 S — 7 s /47rA must be suppressed to zero. Hence, the continuum limit corresponds to the weak 
coupling limit S s — * 0. _ 
From ( 2.39| ), the damping kernel 7(w) in (2.37) is calculated analytically. The leading term is given bytf 

s-l 



7(z) 



7s 



sin(7rs/2) Vcj 



(|) [l + 0{z/^ c ,{z/w c f-% 



7s z 



In 1 



27s 



tt(s - 2) 



(1 + 0(z 2 /col (z/u; c y- 2 )) 



(0 < s < 2), 

(s = 2), 
(s > 2). 



(2.40) 



From ( 2.38| ) and (2.40), it is proved that the delta function disappears at u> — in Re a(u>) for < s < 2. This 
result strongly indicates that the Drude weight D in (2.19) vanishes for < s < 2 at all temperatures and 
damping coefficients. This is because the continuum limit describes the low-temperature and weak-coupling 
limit, T, S s — > 0, while this limit should give the most coherent result in the parameter space. On the other 
hand, the Drude weight exists for s > 2, and the damping effects appears only in the mass renormalization 
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Fig. 2. The optical conductivity <t(uj) in the continuum limit for s = 0.5, 1, 1.5. Particularly, for the ohmic damping 8 = 1, 
fx(oj) takes a simple Drude form. 



In this case, the residual conductivity for uj > vanishes. The above result is related to the Brownian motion 
caused by a heat bath, where the long time behavior of ({q(t) — q(0)) 2 } at finite temperatures is proportional 
to t 2 for s > 2, and t s for s < 2. The transition at s — 2 can be understood by a simple discussion on the 
mass renormalization expressed in general asu'lS* 

M^M+ - P (2.42) 

71" Jo W 

When this integral is convergent, the particle motion becomes coherent in the continuum limit. 

Let us see the s dependence in more detail. The residual part a Tes (u)) appears only for the case s < 2. 
From Q2.3S| ) and Q2.4C| ), we obtain 

\ s-1 

U) 

<Jrcs(u) = —^ — n , (2.43) 

V ' My / , \s-i\ 2 f X 2( S -1)' 

- - cot — I - 

where 7 = (7 s a; 1 ~ s ) 1 ^ 2 *\ In this case, cr rcs (ti;) behaves qualitatively as 

'-(«>-{£: ttl < 2 ' 44 » 

The w-dependence of a res (u) for s = 0.5, 1, 1.5 is shown in Fig. 2. For subohmic damping < s < 1, the 
conductivity is reduced to zero as u> decreases. On the other hand, in the superohmic case 1 < s < 2, the 
conductivity diverges for u> — > 0. For the ohmic damping s = 1, the optical conductivity <7 re s(w) takes a 
simple Drude form. 

As s increases toward 2, the curve of a les (uj) approaches the 1/u) form. The case of s = 2 is marginal, and 
we obtain 

6^72 1 
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The optical conductivity a Ies (u>) behaves nearly as l/u> for to <C u> c . For s > 2, the incoherent part vanishes, 
and the delta function appears in Recr(w). 

2.4 Incoherent tunneling regime 

In the limit of high temperatures and/or strong damping, the particle moves incoherently to the neighbor 
sites. In this regime, the occupation probabilities obey a simple master equation |3lIP and the hopping rate 



is formulated by the Fermi's golden rule, where the i ncohe rent motion of the particle is described by the 
term of order of A 2 in the correlation function. From ( 2.26| )-(2.28), the relevant term is given by 



A(r) = 2A 2 exp [-(f>(r)] . (2.46) 

From (2.30), the real-time correlation function A(t) is derived as 

A(t) = 4A 2 sin R(t) exp(-5(t)), (2.47) 

where S(t) and R(t) are real functions defined by 4>(r = it) = S(t) + iR(t). The explicit form of S(t) and 
R(t) is given from (2.3) and (2.5) as 

dw^^(l - coswi)coth^, (2.48) 



a 2 f c 
k Jo 



R(t) = — I duj^-smujt. (2.49) 



The Fourier transformation of (2.47) gives 



ImA(w)=4A 2 / dt sin ujt sin R(t)exp(-S(t)). (2.50) 
Jo 

Thus, the optical conductivity is calculated from cr res (w) = e 2 o 2 ImA(w) jui. The formulation of the incoherent 
region is valid only when the integral ( 2.50| ) is convergent. 



The above formulation is quite similar to the calculation of the nonlinear mobility in biased periodic- 
potential systemsS Actually, the tunneling rate T(e) and the nonlinear mobility /z(e) with a bias e are 
calculated in these systems as 



/>oo 

r(e)=4A 2 / dt cos et sin R(t)exp(-S(t)), (2.51) 
Jo 

/•OC 

/ dtsmetsmR(t)exp(-S(t)). (2.52) 
Jo 



, s 4a 2 A 2 
H(e) = — — 



Compared (2.52) with ( 2.50| ), the optical conductivity is. obtained as CT ros (w) = e 2 /i(e — > to). Further, it is 



proved that fx(e) is related to r(e) through the relatione^ 

fi(s) = a 2 tanh(/3e/2)r(e) /e. (2.53) 
Thus, we can obtain the optical conductivity from the tunneling rate T(e) as 

arcsH = e 2 a2 tanh(^/2) r(£ ^ 

In order to calculate T(e) in the incoherent regime, we can utilize the results on two-state systems obtained 
in ref. [t], where the tunneling rate is formulated by the same form as (2.51), and is calculated for a few 
regions in the parameter space. For the ohmic damping, it is proved that the expansion in terms of A gives 
the systematic high-temperature expansion. We consider the lowest contribution of order of A 2 in § [3[ In 
the subohmic case (0 < s < 1) and the superohmic case (1 < s < 2), the calculation is so complicated that 



analytical treatment is possible only in the limiting cases. The results are given in Appendix A, and we do 



not discuss the other cases in this paper. We note that the result ( 2.5C ) is also related to the transition rate 
of the photo-induced tunneling in the two- level systems l3 

In the incoherent tunneling regime, the partition function (2.12) is also approximated as 

Z = l + 2/3A 2 / drexp[-0(r)]. (2.55) 
Jo 

This form corresponds to the partition function of dissipative two-iitate systems except a difference of a 
factor 2 which comes from two possible transitions to neighbor sites-O* 

2.5 Weak coupling theory 

For weak damping (5 S *C 1), analytical calculation is possible for all temperatures. We, however, should 
be careful to constitute the weak coupling theory, because we must deal with screening effects of interacting 
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classical particles described in (2.12). In order to treat the screening effects, we adopt the 'ring approxima- 
tion', which is proved to be equivalent to the Debye-Hiickel theory. Since this approximation is believed to 
be valid in the weak coupling region, we expect that reliable results is obtained in the weak damping region 
by this approximation. In this subsection, we derive the partition function and the optical conductivity 
based on the ring approximation. Since details of the calculation are straightforward, but rather tedious, we 
show only the results here. Details of the calculation are given in Appendix [b| and [(]. 
2. 5. 1 Partition function 

First, we define a Fourier transformation of the potential </>(r) as 

rP 

ire iuJmT (j)(T), (2.56) 



where ui m = 2irm/f3 is the Matsubara frequency. By comparing (2.3)-(2.4) with (2.37), (f>(kj m ) can be related 
to the damping kernel j(z) defined in (2.37) for io m ^ as 

<f)(iu> m ) = -—^-%z = uo m ). (2-57) 

To avoid unphysical divergence, we shift the potential 4>{t) by the zero frequency component 0o = </>(i^m = 0) 
as 



<p\T) -> cp{T) -t- cp . 
By the potential shift, the transition amplitude A is renormalized as 



(2.58) 



A = Aexpf-^» ) , (2.59) 

where the factor exp(— </>o/2) corresponds to the Frank-Condon factor. Details of calculation for <fio is given 
in Appendix [b|. 

From the approximate calculation given in Appendix C, the analytical form of the partition function is 
obtained as 

r2-rr 

i i n 

(2.60) 



C c\f) 

Z = J Q 2^ [U{n)] 



U(n) = n+ dnQ(n) -nQ(n), (2.61) 
Jo 

<3W = m 1^ i Tr 2 ' 62 ) 

P ^ n 1 ~ n0(iw m )//3 

where n = n{9) is determined by the equation 

2/?Acos6> = ne~ Q{n) . (2.63) 

From eqs. (2.60)-(^J3S ), the partition function Z is calculated analytically at least for the ohmic damping. 
(See§||) 

2.5.2 Optical conductivity 

The ring approximation is also applicable to the conductivity a(uS). When we define a screened potential 
<p(r; 9) as 

<P(t; 0) = > j^-j. --re m , (2.64) 

lUJrn 

the correlation functions are expressed as 

A(r)=A 1 (r)-A 2 (r) (2.65) 

= h r Hp&e^ ™ + 2ie + ^ 6)] • (2 - 67) 
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After the replacement r — > it, the screened potential is expressed by 

p(ii;0) = S(t ; 0)+iR(f;0), 



(2.68) 



where S(t;0) is a real part, and R(t;9) is an imaginary part. Using this notation, the real-time correlation 
function in (2.30) is written by 

r 2TT 



A(t) 



d6 



M<>) 



2tt (3 2 cos 2 6 



sin R(t; 



(2.69) 



Thus, the optical conductivity is calculated from (2.21) and (2.31). 

In the high-temperature limit (J3 — > oo), the screened potential (2.64) agrees with the unscreened potential 
4>(t), because the second term of the denominator in ( ]2.64 ) is suppressed. Therefore, the expressions of the 
weak coupling region is consistently connected to the formulation in the incoherent regime, ( p. 46 ). 

In the limit T, S s — > 0, the weak coupling must correspond to the continuum limit. To see the relation, it is 
convenient to study the imaginary-time correlation function A(t). In the continuum limit, the dimensionless 
coupling coefficient S s is reduced to zero. Therefore, the screened potential ip(r) becomes small, and the 
approximation 

exp(±v3(r)) ~ 1 ± tp(r) (2.70) 
is justified. From (2.65)-(2.67), the imaginary-time correlation function is calculated approximately as 



A(iw TO ) = - — 



d6 



2n (3 2 cos 2 6 



Me) ^ + e 2w^ 



</>(iw m ) 



n(6)<j)(iLOm)/P 



(2.71) 



The continuum limit (A — > oo) corresponds also to the low-temperature limit (T/A — > 0), where the 
contribution at 9 — is dominant in the integral ( 2.71 ), and n is evaluated by n = 2/3A from ( 2.63| ). Thus, 
we obtain 



A(icj„ 



7(w m ) 



Ma 2 Lo m + j(u) m ) ' 



(2.72) 



where M = 1/2 Aa 2 and the relation ( [2.57 ) is used. From ( 2.72| ), the optical conductivity in the continuum 
limit is obtained as 



a Tes (u>) = e a Im 

= Rc — — 

M z + j(z) 



A(iw m —>■ uj + iS) 



(2.73) 
(2.74) 



This expression corresponds to the optical conductivity in the continuum limit, ( 2.3S ). 

§3. Ohmic Damping 

3. 1 General review of ohmic damping 

In this section, we consider the ohmic damping case, in which the spectral density is described as 

2nK 



-cue 



■j/u 



(3.1) 



Here, if is a dimensionless coupling constant. The Fourier component of the potential in (2.3) is calculated 

as 

2nK 



<Xiw m ) 



2JHog 



2tt 



(0 < \w m \ < u c ), 
(u> m = 0). 



(3.2) 



The calculation of (f>Q = <j)(iu) m = 0) is given in Appendix [b|. From ( |3 . 2| ) , the renormalized hopping amplitude 
A in ( 2. 59] ) is given by 



(3u c 
2tt 



—K 



(3.3) 
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Ik Incoherent y Solvable line K=1/2 
i\ ^ (Toulouse limit) 

V 



Incoherent 




localized at T=0 



Fig. 3. The phase diagram of dissipative two-state systems for the ohmic damping. This phase diagram has been obtained by 
the Non-Interacting Blip Approximation adopted in rcf. ^. In the 'coherent' region, the expectation value of the position, (q(t)) 
shows oscillations, while in the 'incoherent' region, (q(t)) shows the exponential decay. The regions which allow analytical 
treatments is also shown: the weak coupling regime and the solvable line (K = 1/2). The localization-delocalization transition 
at T = is also drawn. 



In the following subsections, we calculate the specific heat and the optical conductivity for the ohmic dissipa- 
tion. Before showing the results, we briefly review previous works on the properties of the ohmic dissipation. 

§or two-state systems, nonequilibrium behaviors of the dissipative particle have been studied by Leggett et. 
They have calculated the real-time evolution of the expectation of the particle position with the Non- 
Interacting Blip Approximation (NIB A). The properties of the system are determined by the damping 
strength K and the temperature T, where the phase diagram obtained by NIBA is given in Fig. 3. There 
exist two regions in the phase diagram. In the 'incoherent' region, the particle moves incoherently, and 
the expectation value of the position, (q(t)} decays exponentially to the thermal equilibrium state. On the 
other hand, in the 'coherent' region, (q(t)) shows damped oscillation. We should note that several different 
definitions of the terms, 'coherent' and 'incoherent', arfiLpossibleJlS' For example, in the numerical study of 
equilibrium properties of dissipative two-state systems,EfEl' the term 'incoherence' is defined by disappearance 
of an inelastic peak of the response function, which is observed at K > 0.33. A different definition of the 



term, 'incoherence' is also used in § 3.4 



As compared with the two-state systems, it is more difficult to solve the lattice model studied in this 
paper. For example, the simple approximation such as NIBA is not applicable to the lattice model. The 
dissipatiy^pacticle in the tight-binding model has been studied within the real-time path integral formulation 
in detail£3'EirE-$ Particularly, early discussions are based on the the duality transformai | ipa |r attliich relates 
the tight- binding model to the periodic-potential system in a weak corrugation limit .E3cSe3^E3) In these 
works, only the mobility of the dissipative particle has been of interest, and the whole dynamical properties 
of equilibrium states such as the optical conductivity have not been considered. Further, the mobility at 
intermediate temperatures has not been obtained exactly except for the special value of the damping strength 
(K = l/2). 

In this paper, we consider mainly two limiting regions in the phase diagram. The first one is the weak 
damping region K -C 1, which allows analytical treatment at all temperatures. The formulation of the weak 
damping theory can connect the low-temperature 'coherent' region to the high-temperature 'incoherent' 
region in a natural way. The result obtained in this narrow region will also be useful to understand properties 
in the other region. The second region is the high-temperature and/or strong damping region. In this limit, 
the dissipative particle moves incoherently with the hopping rate determined mainly by the Fermi's golden 
rule. Especially, in the region K > 1, the high temperature expansion is convergent, and the Fermi's golden 
rule gives the leading term at all temperatures. On the other hand, for K < 1, the high temperature 
expansion fails below the Kondo temperature Tk- The summary of the results on these two region is given 
in Table I. 
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Low-temperature and 
weak-damping region 



High-temperature 
or /and strong damping 



Specific heat C 


oc T/K 


H 


r r- 


-2+2X 


fif < 3/2) 
(K > 3/2) 


Optical conductivity tr(oj) 


x f u>-*+™ (^» 70 ) 
\ const. (u < 70) 


-] 


i T " 


2+2K 
-2+2K 


(w > KT) 
(w « iCT) 


DC conductivity ctdc 


oc 1 — const. T 2 






ocT- 2 - 


-2A" 



Table T, Summary of the results. Here, K" is the dimensionless damping coefficient, and 70 is the damping frequency defined 

by ( 3.34 ) . The specific heat at high temperatures is obtained from the calculation in two-state systems in ref . [39|. 



With regard to the optical conductivity, the result on the continuum limit is also useful to understand 
the general behavior. The continuum limit corresponds to the low-temperature and weak-coupling limit 
(K,T — ► 0), and the particle is expected to move coherently. From ( 2.43 ), the optical conductivity in the 
continuum limit is calculated as a simple Drude form 
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e 7 



M{lu 2 +7 2 )' 



When the finite damping is introduced, the optical conductivity deviates from the Drude form (| 
behavior is studied in detail in S 3.2. 



(3.4) 



4|). This 



In addition to the above regions, it is expected that the calculation at K = 1/2 is tractable analytically 
by following ref. though we do not deal with the case in this paper. In the other regions (particularly 
for < K < 1), the analytical method cannot be used. Hence, the dynamical properties of these region will 
have to be studied by numerical calculation, which has not been performed before to the lattice model to 
our knowledge. 

Finally, we comment on the renormalization of A due to the damping. For the ohmic damping, the relevant 
frequency of the system is given by 

K/(l-K) 



/AY 



(3.5) 



for < K < 1. The frequency gives the scale of the Kondo temperature Xk- By rewriting A e g- as A, the 
final results can be expressed without the cutoff frequency w c .0' In this paper, all the results are expressed 
by A off finally. 



3.2 Specific heat 

In this subsection, we consider the specific heat of th e sys tem coupled to the ohmic heat bath. At high 
temperatures, we can obtain the specific heat C from ( 2.55| ). This result is exactly twice as large as the 
one obtained in dissipative two-state systems. As a result, the specific heat is proportional to T 2K ~ 2 for 
< K < 3/2, and proportional to T for K > 3/2. For details of the cal culat ion and result, see ref. |3^. 

Next, we focus on the weak damping region K <C 1. From eqs.(2.60)-( 2.63 ), we obtain 



Z = 



2 " M u {t 
2tt 



U{n) =n + log T{Kn + 1) - Kmp(Kn 
Q(n)=K(ip(Kn + l)+*f), 
2/3A cos<9 = ne- Q[n \ 



1), 



(3.6) 

(3.7) 
(3.8) 
(3-9) 



where 7 is the Euler's constant, and ip(z) is the polygamma function. From this expression, the specific heat 
is obtained as Fig. 4 by numerical integration. We can see that, even for the weak damping, the specific heat 
is suppressed at low temperatures, and goes to zero in the limit T — ► 0. This behavior is different from the 
dissipationless case K = 0, where C approaches a nonzero constant value at low temperatures. This result 
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Fig. 4. The specific heat for weak damping (K = 0.03,0.06,0.1,0.15). The specific heat is suppressed to zero as the temper- 
ature decreases. This behavior is quite different from the dissipationless case (K = 0), where the specific heat approaches a 
nonzero constant at iow temperatures. 



can be interpreted as follows. In the dissipationless case, the partition function is classical in the momentum 
in the sense that the Hamiltonian is diagonalized in the momentum space. We note that the quantization 
of the momentum due to a boundary condition is not taken in this formalism so that the specific heat 
approaches a classical nonzero value. When damping is introduced, the Hamiltonian cannot be diagonalized 
in the momentum space, and the quantum effects suppress the specific heat at low temperatures. 
The asymptotic expansion of Z for (3 — 1/T(^> 1) gives 

In Z = 2(1 - K)PcA cS + \ InK + + {P ~% (3.10) 

where c is a renormalization factor defined by c = (4:7rK^) K ^ 1 ~ K \ and takes c ~ 1 for weak damping 



A' « 1. Details of the derivation is given in Appendix D. From the partition function, we obtain the 



low-temperature behavior of the system energy E and the specific heat C as 

E = - A On*) = -2(1 - K)cA cS + + OCr .) (3 .n) 

pi rp rp 

We can see that the specific heat is proportional to the temperature T, and the coefficient depends on the 
damping strength as C ~ T/KA c g. We denote the relevant energy scale with 7 = 4irKcA c g. As shown 
later, 7 expresses the relevant energy scale in the optical conductivity for weak coupling region. Then, we 
can write the specific heat as 

C ~ 7rT/3 7 (3.13) 

at low temperatures. 

Next, we discuss the density of states, D(u>) of the dissipative particle. From the partition function, D{u>) 
is calculated by 

/oo 
AioD{Lo)e~^. (3.14) 

Here, the energy is shifted in order that the ground state energy becomes zero, and the infinitesimal quantity 
8 is introduced for convenience. In the dissipationless case, the density of states of ID lattice systems is 
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Fig. 5. The density of states in the dissipative case (thick line) and the free case (thin line). By introducing dissipation, it 
seems that the states expressed by the patched area moves to the weight of the delta function at u) = 0. 



calculated as 



D(u) = 




(3.15) 



W 

Y 



where W is a band width. Particularly, for low-energy states, D(cu) behaves as 



-1/2 



(3.16) 



In the presence of dissi patio n, the density of low-energy states is strongly modified. To see this, we rewrite 
the partition function ( |3 .10 ) by using the effective band width W = 2cA c g as 



K e l/12KW 

Vk( i 



1 



12KW 



0(/3- 



Then, the low-energy form of D(u>) is obtained from ( p,14| ) as 

1 



D{oj) ~ VK6(u) + 



(3.17) 



(3.18) 



The difference between ( 3.16 ) and ( 3.1§| ) is schematically shown in Fig. 5. We can see that the dissipative 
environment strongly modifies the low-energy states. Roughly speaking, because of dissipation, some parts 
of low-energy states up to uj ~ KW (the patched region in Fig. 5) gathers to to — in the form of the 
delta function. Although we cannot discuss the details of this change until the one-particle Green's function 
is obtained, it can be conjectured that the dispersion is modified to be flat by dissipation, and the wave 
function of the low-energy states are localized in space. 

We note that similar behaviors have been shown in two-state systems, .-where the Shottkey form of the 
specific heat is modified to the T-linear behavior at the low temper atures.ES Also in this case, the T-linear 
behavior of the specific heat is induced by the modification of density of states caused by the environment. 

3.3 Optical conductivity 

3. 3. 1 The limit of high temperatures and/ or strong dissipation 

For the ohmic damping, S(t) and R(t) in (2.48)-(2.49) are obtained as 



S(t) = 2K\og ^ 
R{t) = ?rifsign(t). 



(3.19) 
(3.20) 
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Fig. 6. The optical conductivity cr res (w) at high temperatures (T = 6A e ff , 12A e ff , 18A e fj) for K = 0.5. The dashed gray line 
represents the asymptotic behavior (oc l/io) at u>/A c g S> 1. The characteristic energy at which the optical conductivity 
deviates from the asymptotic form is estimated as KT. 



From fl2.50| ), we obtain the optical conductivity a les (uj) = e 2 a 2 ImA(w) / 'u> as 



" J — I P u c 
CTres(w) = 2e a — — — 



1-2K 



2e J crA cff 



uj c \ 2ir 
fPA cS 



V 2tt 



\T(K + i/3w/27r)| 2 sinh/3cj/2 
T{2K) u 

K \T{K + i(3uj/2tt)\ 2 sinh/?w/2 
T(2K) u> ' 



(3.21) 



(3.22) 



The temperature and frequency dependences of er rcs (w) are shown in Fig. 6. In the low- frequency side 
(uj -C KT), a Tes (uj) depends on the temperatures. The DC conductivity aoc = c rcs (a; — > 0) is calculated as 



gdc = 2ne 2 a 2 



f(3A cS 



2-2K 



my 

T(2K)'- 



(3.23) 



V 2tt 

which is proportional to T 2K ~ 2 . On the other hand, in the high-frequency limit, a les (uj) takes a temperature- 
independent form 

2-2K 



2-Ke 2 a 2 ( A. 



T(2K) V w 



Vir 



(3.24) 



On the high-frequency side, a ICS (w) is proportional to oj~ 2+2K , and deviates from the Drude form (3.4). 
3.3.2 Weak coupling region 

Next, we consider the weak damping region, K « 1. For the ohmic damping, the screened potential 
defined by (2.64) is calculated as 



<p(t\9) = 



-2-kK 



where 7(f) is an inverse of a screening length in the r-direction given by 

7(0) = 2nKn{9)/(3. 



(3.25) 



(3.26) 



Following the usual way, the sum over the Matsubara frequency in (3.25) is replaced by an integral form. 
As a result, ip(it; 9) — S(t; 9) + iR(t; 9) is obtained as 



R{t;6)=irKe-^ t 



(3.27) 



S(t; 9) = 2K 



dw 



u 2 +l(9f 



2 Buj 
coth — cos 

(3lu 2 
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Fig. 7. The optical conductivity cr(w) at zero temperature for weak damping (K = 0,0.1,0.2). As the damping increases, 
cr(w) deviates from the Drude form (K = 0). At high frequency, cr(u>) is proportional to u> 2K ~ 2 . 



2 , Plu 

— coth — — COS U)t 

(3lj 2 



(3.28) 



where 0(x) is a step function. For no damping case 7 -* 0, the above expressions become unscreened 
potentials, S(t) and R(t) in (2.48)-(2.49). Analytical expressions for S(t;9) are obtained only for two 
limiting cases. At high temperatures (3 — » 0, S(t; 9) is calculated as 



On the other hand, at low temperatures, S(t; 9) is expanded by the temperature T = 1//3 as 



S(t;0) = S (t;6) 



C(T 2 ), 



S (*; A) - K \e^ 0)t Y l i{~ 1 {9)t) + e-^ t Ei( 1 (9)t) 



(3.29) 



(3.30) 
(3.31) 



where Ei(z) is the exponential integral function, and So(t) is a temperature-independent part. In both cases, 
the limiting value at t — > 00 is given by 

„s 27ri4:T 



5(i 



cxd; 1 



1(0) 



(3.32) 



Although we can calculate a(cj) by the expressions (3.27)-(3.28) for all temperatures, we only consider 
zero temperature in this paper. At zero temperature, the dominant contribution comes from 9 = in the 



integral (2.69). The real-time correlation function for K <C 1 is obtained by using the asymptotic form 
n = 2/3cA e ff as 



AM 



7o 



die^-^cosh [S (t)]. 



2izK j 

Here, 70 = j(9 = 0) is the damping frequency at T — given by 

70 = AnKcAes, 



(3.33) 



(3.34) 



where c = (4TrKe^) K ^ 1 ^ K \ The optical conductivity can be calculated by numerical Fourier transformation 
of ( 3.33] ). The result for cr rcs (a;) is shown in Fig. 7. The Drude formula (3.4) in the continuum limit is also 
shown. At the high frequency side, the optical conductivity of the dissipative system deviates from the 
Drude form. 
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The qualitative behavior of <j{u>) is obtained by considering the limiting from of So(t): 



So(t) 



2K\og(e^ t) (t 
2K 



0), 



(3.35) 



where 7 is the Euler's constant. From ( p.33| ) and ( p.35| ), the limiting behavior of the optical conductivity is 
obtained as 

r, 2 2 / a \ 2-2K 

, [oj > oj*), 



a(oj) 



T{2K) 



To 2 



(3.36) 



2irK 1I+UJ 21 



(oj < OJ*), 



where oj* = (2c) 1 ' 2K is the crossover frequency at which the two expressions take t he s ame value. The 
high-frequency behavior of a(oj) agrees with the result of the high temperature limit ( 3.24 ). This i ndicates 
that the high-frequency side is independent of the temperature at any damping strength, since ( 3.24 ) is valid 
at any value of K. On the other hand, the low- frequency side coincide with the Drude form (3^) obtained 
in the continuum limit. It is interesting that the asymptotic form in the limit oj — > 00 of the Drude form 



(the first equation in (3.36 ) ) d eviates from the high frequency form (the second equation in (3.36) ) shown 
in the lower equation of (3.36) by factor 2. At first glance, one might have suspicion on this difference. 
However, this difference indeed exists: When we take the limit K — > 0, the crossover freq uency oj* in ( 3.36 ) 
goes to infinity. Then, the region w > w* disappears and the lower expression in ( |3.36| ) governs even the 
high-frequency region, which reproduces (|3.4|) . 

On the high frequency side w > w*, the particle moves by using the heat-bath excitations. Therefore, 
the particle, moves incoherently even at zero temperature. The results are relevant to the photo-assisted 
tunneling .liBP 
3.3.3 DC conductivity 

We consider the DC conductivity ctdc of the lattice model for the weak coupling region (K <C 1). The 
low-temperature behavior is obtained by the asymptotic expansion over (3=1 /T as 





DC 



'"Kss) + °< T3 >' < 3 - 37 > 

e 2 a? /2wK. (See Appendix Q|.) Thus, the thermal fluctuation suppresses ctdc- Tb£_T 2 - 
suppression at low temperatures is consistent to the general study of the linear mobility fii = ctdc / e 2 £5 Jn 
the previous work, the linear mobility has not been obtained correctly even in the weak damping K <C 1 .lj 
The high-temperature limit of ctdc for K <C 1 is obtained from (3.23) as 



where CTq C = 



cdc/odc = 8-n- 2 



gAeg 

2vr 



2-2K 



(3.38) 



We compare this result with the linear mobility hi obtained in ref. 33. For the special value K = 1/2, the 
linear mobility is calculated for all temperatures as 



2T 



-V>'(A cff /2T + l/2), 



(3.39) 



where ip'(z) is the trigamma function, and /io = a 2 /2irK. At low temperatures, the DC conductivity 
°dc = /"(/e 2 behaves as 

T 2 



o'dc/o'd 



DC 



1 



3A? 



(3.40) 



and show the T 2 -suppression. This strongly indicates that the T 2 -suppression appears for any value of 
< K < 1/2. At high temperatures, the DC conductivity for K = 1/2 shows a 1/T behavior as 

a DC /a« c ~ (3.41) 



which agrees with ( 3.23 ) for K = 1/2 
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3-4 Application of the results to real systems 

The tight-binding model of a particle with ohmic damping is relevant to following systems: the muon 
motion in a metal, the incoherent carrier motion in coupled Tomonaga-Luttinger liquid and Josephson 
junction systems. We, however, must comment on the qualification in applying the result obtained in this 
section to each system. 

The mu an .motion in a metal has been discussed first in connection with the Anderson's orthogonality 
theorem,oE§E!J and subsequently this problem has been formulated based on the path integral.Ef Exper- 
iments in Cu have showed the T 2K _1 -behavior in the diffusion rate as expected in the incoherent regime 
of the Caldeira-Leggett model. In this experiment, the effective hopping amplitude A c «_~ O.lmK is so 
small that the low-temperature region is difficult to be observed in the presence of disorder S3 Moreover, the 
Caldeira-Leggett (CL) model is thought to fail in describing the heavy .particle motion at low temperatures, 
since the damping is overestimated as pointed out by several authors-C^'Ed' The authors have showed that 
the damping strength due to the conduction electrons becomes small at long distances, while in the CL 
model, this feature is absent, and the damping is overestimated at long distances. This difference would 
be critical particularly at low temperatures, since the configuration paths with large spatial distribution are 
significant. 

Recently, the CL model has been used to discuss the incoherent non-metallic transport in strongly corre- 
lated anisotropic metalo-* Based on two coupled Tomonaga-Luttinger(TL) liquids model, the authors have 
discussed inter liquid particle motion with small hopping matrix clement, by taking a dissipative two-level 
system as an effective model. It is, however, problematic to extend this treatment if infinite number of TL- 
liquids are coupled. The effective Hamiltonian does not correspond to the Caldeira-Leggett model, because 
the damping coefficient again should depend on distance between chains. Hence, it seems also problematic 
to take the naive CL model as an effective model of coupled TL-liquids for calculating low-energy properties 
at low temperatures, and the application of the tight-binding model for incoherent interchain hopping should 
be considered more carefully. _ _ _ 

The dissipative phase motion in Josephson junctionsE-i E3'E3* seem to be free from the problems in contrast 
to the previous two systems, and the result obtained in this section may be applicable. In this system, 
however, damping is so strong (K ~ 1) that the coherent motion has not been confirmed yet._In this 
system, the incoherent tunneling, motion with the rate proportional to T 2K ~ 1 has been observed,© while 
macroscopic quantum coherenceLP in the two-state systems expected in the weak-damping region (K <C 1) 
has not been observed so far. 

Thus, unfortunately, no system is at the moment available to directly apply the result in this paper. We 
wish that progresses in experimental techniques would make it possible to realize experiments to which the 
result in this paper can be applied. 

We also speculate that the present model may be applicable to the incoherent carrier motion in strongly 
correlated electron systems. For example, the retracable path approximation for correlated electron sys- 
temsBEj'E-P gives the optical conductivity of non-Drude forms roughly proportional to 1 jui in contrast to 
the Drude form a(u>) ~ l/(^ 2 + 7 2 )- This behavior can be reproduced in the Caldeira-Leggett model for the 
ohmic damping case (s = 1,K ~ 1/2) or the superohmic damping case (s ~ 2). In the retracable approxima- 
tion, the coherent motion of the particle is strongly suppressed by the antiferromagnetic background. This 
mechanism is similar to the present model of the dissipative particle, when the heat bath is regarded as the 
antiferromagnetic background. The similarity is, however, not complete, since the density of states, D(u>) 



for the ohmic case given in (3.18) differs from the result obtained by the retractable path approximation 
at low-energy states. Nevertheless, we speculate that the lattice model of the dissipative particle may be 
utilized as a basic model of the decoherence of the particle due to the electron-electron correlation at least 
in the high energy region of the optical conductivity. 

§4. Summary 

We have studied the thermodynamics and transport properties of the dissipative particle in the tight- 
binding model. By utilizing the imaginary-time path integral, the specific heat and the optical conductivity 
have been formulated for arbitrary form of the spectral density. A systematic approximation has been 
considered to treat the weak-coupling region at all temperatures. We have obtained an analytical form 
which can connect the high-temperature region to the low-temperature region. 

The actual calculation has been performed for the ohmic damping case. The specific heat at weak damping 
shows T-linear behaviors at low temperatures, and the density of low-energy states is modified by dissipation. 
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The optical conductivity shows non-Drude from even at zero temperature, and behaves as lo 2K ~ 2 on the high- 
frequency side, where K is the dimensionless damping strength. At high temperatures, the high-frequency 
side does not change the non-Drude form, and only the low-frequency side depends on the temperature T. 
Particularly, the DC conductivity is proportional to T 2K ~ 2 . We also comment on the application of the 
Caldeira-Leggett theory to real systems including the descriptions of the incoherent transport in correlated 
electron systems. 

In this paper, we have focused only on the ohmic damping. The detailed calculation for the non-ohmic 
damping and the solvable line K — 1/2 remains for future study. We expect that the simple model considered 
in this paper may contribute to understand transport properties of a particle coupled to other degrees of 
freedom, and wish that the results of this paper may be directly observed in a properly fabricated experiment. 

Appendix A: Non-ohmic Damping in the Incoherent Tunneling Regime 

In this appendix, we study the optical conductivity <r(w) for the non-ohmic damping case in the incoherent 
regime following the formulation in § 2.4. The results on the tunneling rate T(e) obtained in ref. [?| are utilized 
to obtain a(ui). We only consider a few limiting cases which allow analytical treatment. 

We begin with the superohmic damping 1 < s < 2. Under the strong bias e ^> A and the weak coupling 
condition 5 s {e/ui) s ~ 1 <C 1, the tunneling rate T(e) is calculated as 

a 2 /2A\ 2 

r(e)-yf— J J(£)coth(/^/2). (A-l) 



By using (1.8) and (2.54), we obtain the optical conductivity as 

2 

/ \ \ 

cr(w) 



(A-2) 



which is temperature-independent at all temp eratu res. Hence, the optical conductivity behaves as a(u)) cx 



. This behavior corresponds to the result ( 2.44 ) in the continuum limit (T, S s —> 0) on the high frequency 



side, except a difference by a factor 2. Similar difference appears also in the case of the ohmic damping, and 
the origin of the factor 2 is expected to be the same. (See ( 3.36 ) and the following text.) 
The tunneling rate for e = is calculated as follows: 



(2- s) sin(7rs/2) 



1/(2-*) 



3- s 



A2 

UJ 



1/(2-*) 



T* = 



(A-3) 



(A-4) 



The DC conductivity is obtained from ( |A-3| ) by odc = e 2 a 2 T/2T, and behaves as cjdc oc T~^ 7 at high 
temperatures. The exponent of T decreases from —2 to — oo as increasing s from 1 to 2. 

Next, we consider the subohmic damping < s < 1. At zero temperature, the tunneling rate is calculated 
under the condition S s (e/uj) s ^ 1 ^> 1 as 

1/2 



2S s T(s - l)sin7r(s - 1) J \2 - s 

This result is valid at high temperatures T* <C T(<C w c ), where T* is given by 

A /£ 

<J s r(a-i) \K 



T(e) = 



(BA) 2 f2ir[26 s T{s)} 



1/2 



4a; 



(l + s)/2s 



x exp 



1 -s 



\2S s r( s )} 
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(l-s)/s 



(A-5) 



where B = cxp[S s \T(s — l)|(cj/a; c ) 1_s ]. This result cannot be related to the result of the continuum 
limit because of the strong coupling condition S s (e/uj) s ^ 1 ^> 1. The optical conductivity is obtained by 
<j(uj) — a 2 e 2 r (w) / 'uj , and behaves as <j{oS) ~ uj^^ 1+3s ^ 2s exp(— const. a;^ 1- ^/ 5 ). Notice that the optical 
conductivity vanishes with an essential singularity as e — > 0. 

The tunneling rate for e — for the condition 5 s (ui /T) 1_s ^> 1 is given by 

(l+s)/2 / \ 1/2 

i i ) \ i / . i . \ 
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2(1 + s) 2 r( s )C(l + s)S s 
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From (A-6), the DC conductivity is obtained by ctdc = a 2 e 2 T/2T, and behaves as <tdc ~ 
j 1 -(3+s)/2 exp(— const, x T s_1 ). The DC conductivity vanishes with an essential singularity as T — ► 0, 
and this result is consistent with the continuum limit, where <7 rog (cL> — > 0) = 0. Note that the power law 
cdc ~ T 2K ~ 2 for the ohmic case s = 1 is obtained from (A-6) by carefully considering the limit s — > 1~. 

Appendix B: Calculation of 0o 

In this appendix, we evaluate the zero frequency component 
4>o is explicitly given as 

„2 



0). From (2.56) and (2.3)-(2.4), 



'jo = 



= — da; 



J(u) 



coth — 

2 /3uj 



We consider the analytical form of the spectral density, 

2ttS s 



J(cj) 
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wher e the exp onential cut-off is adopted, and the cutoff frequency uj c is taken as A, 1/(3 <C w c . By substituting 
(jgjg) to pl|), we obtain 



= 25. 



a*.(f)'r W 
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where z = s — 1, and C( z jfl) is the generalized zeta function, and £(z) is the zeta function. In the second 
equation, we have left only the leading terms under the condition (3lo c 3> 1. For superohmic case z > 0, the 
first term in the bracket is dominant, and <po is independent of temperatures. Further, the result is expressed 
generally reproduced by the Frank-Condon factor calculated as 

7T ,/n UJ Z 



when the integral is convergent as in the case of the superohmic damping s > 1. On the other hand, for 
subohmic damping z = s — 1 < 0, the second term in the bracket in (B-4) is dominant. Hence, </>o depends 
on temperatures. For the ohmic case s = 1, the expression (B-4) fails, and careful treatment of (B-3) is 
required. As a result, the leading term in <j)Q depends logarithmically on temperatures as 

(3uj c 



2K log 



2vr 



where K = 8\. In this case, the renormalized matrix element ( [2.59 ) is given by 



A = A 
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(B-7) 



At a first glance, one might regard that this result is unphysical, because the matrix element is reduced to 
zero as the temperature decreases even for the weak coupling region if« 1. We, however, show that the 
renormalization factor (f3ui c /2ir)~ K is canceled by other factors for observables in the weak coupling theory. 

Appendix C: Ring Approximation 

In this appendix, we give details of the calculation in the weak coupling theory (§ for the arbitrary 
form_pf J(u>). We first study the partition function based on the cluster expansion of classical imperfect 
gas.E3) The partition function and the optical conductivity are formulate d by the ring approximation. 

We begin with the partition function (2.12). After the potential shift ( 2.58 ), we obtain 
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x exp 



{<l>(n - r fc ) + (j)(pi - p k )} -EE 0( Tl ~ P k ^ 



(C-1) 



Lk<l 



1=1 fc=l 



where A is a renormalized transition amplitude defined by ( 2.59| ). Here, new integral variables, {77} and 
{pi} describe the positions of positive (£; = +1) and negative (£; = —1) charges. Note that the partition 
function is multiplied by the factor 2m\/m\m\. This factor comes from the number of the ways in which 
2m ch arges are divided into two groups consisting of m charges to guarantee the electroneutrality condition 
( 2.13 ). Then, we consider the cluster expansion for the partition function (CT). To simplify the expansion, 
we rewrite it as the following form 



Z =E 
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Here, the potential 4>{t) is replaced by Mayer's functions 

f++ — 4>( T k-n) _ 1 
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denotes the average 
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Each term in the expansion on fki in (C-2) is expressed by a product of integrals of clusters. We introduce 
cluster integrals 



bt ' k ~M E (II ^T^ 



(C-7) 



clusters i,j 



where ^clusters denotes the sum of all clusters including k positive charges and I negative charges. The 
products are taken over all bonds combining z-th and j-th charges in each cluster, where <Ji and <jj are signs 
of charges. For example, the explicit form of bi t k for small k, I is given by 
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and so on. The above integrals can be viewed by graphical representations as shown in Fig. 8. By the cluster 
integrals bi t k, the partition function is expressed as 
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where mi t k is the number of clusters including I positive charges and k negative charges. The sums and 
products are taken over I and k except for the case (I, k) = (0,0). The prime in {m^}' denotes the 
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b2,o4 



b3,0 = 

b2,i = 



bi,i = 



12 12 

AAA 

12 12 13 



Fig. 8. The diagrams representing the integrals bi m in (C-8)-(CT2). The closed (open) circles denotes the positive (negative) 
charges. Note that 6; m = fe m ,;. 



summation in accordance with the conditions 

^/mi, fe = m, (C-14) 



l.k 



^ fcm;^ = m. (C-15) 



The first (second) factor in the sum J2{ mt k }' m (^'13) gives the number of ways in which m positive 
(negative) charges can be distributed into the clusters. The third part is the relevant cluster integrals, and 
the forth part is the number of ways to combine m^k groups of positive charges to m^. groups of negative 
charges. 

From (C-14)-(CT5), we obtain the constraint on the total number of charges as 

^2(1 + k)m l>k = 2m, (C-16) 

l,k 

and this constraint can be removed by the sum over m. On the other hand, the electroneutrality condition 

5^(1 - k)mi, k = 0, (C-17) 

Lk 



remains. To remove the constraint ( CT7 ), we introduce an integral over a new variable 9 as 

(CMS) 

From (CT3) and ( |CT8 ), the partition function Z is obtained as 

(C-19) 



E E (•••)= E /"f exp(i^(/- fc)m J (-..). 

m=0{m,, fc }' {mi, k } J0 \ l,k j 



For the dissipationless case (K = 0, A = A), bi t k vanishes for I + k > 2 and the partition function is exactly 
calculated as 

Z= — exp(2/?Acos(9) = I {2f3A), (C-20) 
Jo 2lT 

where Io(z) is the Modified Bessel function. Of course, this result is also obtained by usual treatment for 
the tight-binding model of a free particle. 

So far, the obtained expression for Z is exact. In the following discussion, the so-called 'ring approximation' 
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(a) 




(b) 



V 




(c) 



Fig. 9. The diagrammatic graphs representing the expansion in (C-21) and (C-22) are drawn in (a). The graphs in both (b) 
and (c) are representative cluster integrals: (b) each charge has even bonds, while (c) includes charges with odd bonds. 



is introduced. First, we expand L * 3 as 



r + = r 
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(C-21) 
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Graphical representations of the expansion is given in Fig. 9(a), where the bond denoted with n thin lines 
denotes i>(ji — Tj) n . Here, we assume that (I) the dominant contribution comes from the clusters in which 
each charge is combined by 'even' bonds to other charges. Based on the assumption (I), we neglect other 
clusters which do not satisfy this condition. For example, we consider the graphs as shown in Fig. 9(b), and 
neglect the graphs as shown in Fig. 9(c). The assumption (I) is justified in the ring approximation denoted 
later in detail. At this stage, however, in order to avoid confusion, we do not justify the assumption (I). We 
mention this point after introducing the ring approximation. 

The diagrams as shown in Fig. 9(b) have a crucial property that the cluster integral on the graphs remains 
constant when the sign of each charge is reversed. As a result, the cluster integral l\k\b\ t k depends only on 
the total charge number I + k. In this situation, l\k\bi t k c & n b e den oted with (I + k)\bi+k, and the sum over 
I, k under the constraint I + k = m in the partition function ( C-19| ) gives 
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Here £ = 2/3 A cos 8 is regarded as a fugasity of 'one-component' classical charges. Because the remaining 
calculation follows the usual procedure of the classical imperfect gas as given in ref. |5(| we briefly note the 
results. The diagrams, in which all charges are more than singly connected, are called irreducible diagrams. 
The irreducible integrals (3i are defined on the irreducible diagrams with the size {I + 1) as 
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where the product is taken over all bonds in a cluster, and riij is the number of bonds combining i-th and 
j-th charges. Every cluster integral 6; can be expressed as a sum of terms, each of which is a numerical 
coefficient multiplied by a product of powers of the reduced integrals [3i as 

^3 = \0\ + \02 
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In general, it is proved that the equation for bi is 
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where the prime in {m/c}' denotes the summation subject to the constraint 
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We also define a particle density n by 



= ^2 mb m C 
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This value corresponds to the average of the hopping number in the imaginary-time path q{r) in the original 
action, and qualitatively gives an effective hopping amplitude through ~ f3/n. By using j3i and n, the 
partition function Z can be expressed in a closed form. To see this, we expand £ by n as 



£ = a\n + a 2 rL + ■ 



and determine a\, 02, • • • so that (C-28) is satisfied. Thus, we obtain 

a,i = b^ 1 = 1, 03 = 8&2 — 3&3, 

a 2 = -2& 2 , a 4 = -406i + 30b 2 b 3 - 46 4 . 



From ( C-25 ), these results can be rewritten in terms of the /3/'s as 

a-i = 1, 03 = -(/?2 - a-/ 3 ?): 

« 2 = -/?i, a 4 = -(/3 3 -/3i/32 + |/3 1 ). 

In general, £ is expressed by the particle density n as 
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I t can be easily checked that the first four terms in ex pansi on of (|C-32|) gives ( C-31 ). The partition function 
( C-23| ) is expressed in terms of /3; and n by utilizing ( |C-29 ) and ( |C-3l| ). As a result we obtain 
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At this stage, we introduce the approximation by assuming that only the ring irreducible diagrams as 
shown in Fig. 10 are dominant. This approximation is called the 'ring approximation'. Then, the irreducible 
integrals (3i is expressed by the Fourier component <f){iuj m ) as 



2 = 1 



the relation (C-32) and the exponent U(n) is expressed in simple forms as 

£ = ne- Q{n) 

U = n+ [ dnQ- nQ. 



(C-38) 
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Fig. 10. The ring irreducible graphs considered in the ring approximation. 



and thus we obtain 

Q -pj^l-r*<M/P' (C39) 

Here, the sum is taken over the Matsubara frequency, oj m = 2nm/ (3. The justification of the ring approxi- 
mation is not clear. It, however, can be proved that this approximation is equivalent to the DebiiE-Huckel 
theory, which has been considered to calculate thermodynamic quantities in strong electrolytes. lZP In the 
Debye-Hiickel theory, the screening effects due to the other c harg es are described by a screening potential 



ip(r). (The explicit form of the screening potential is given in ( |2.64| ).) These screening effects are essential to 
study thermodynamic quantities at low-temperatures in the original model. In the ring approximation, the 
assumption (I) adopted in the previous discussion is automatically satisfied, because the neglected diagrams 

by the assumption (I) are also removed in the ring approximation. 

We study the condition, in which the ring approximation is justified. From ( C-39| ), it is expected that 



the ring approximation is certainly valid as long as the expansion quantity n<f>(u} m )/{3 is small enough at 
i^m ~ 1/"?) where f ~ /3/n ~ A~g is a average distance between two particles. For the exponent < s < 2 
in the spectral density, this condition corresponds to 

Sa (^r)' 1<<L (c ' 40) 

In the ohmic damping case (s = 1, 8 S — K), this condition is nothing but the weak-damping condition 
K -C 1, where we expect that the approximation gives reliable results in the weak-damping region. 

Next, we calculate the optical conductivity. The correlation function Ai(r) in (2.27) is expressed by the 
cluster integrals as 



Ax (r) = I " g exp ( £ 6 a (/3A)<+<V*«- 

^B lk {f3h) l+k e ie{l ~ k] . (C-41) 
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l,k 

Here, Bi^ is the cluster integral of the diagram including two 'fixed' charges at r' = 0, r, and I and k are the 
numbers of additional positive and negative charges. Note that diagrams in which the two 'fixed' charges are 
separated into two clusters do not contribute to the optical conductivity, because such terms are independent 
of T. 

In the ring approximation, (C-41) is reduced to the partition function of one-component charged particles. 
Thus, we obtain 

Aiw - I ^ exp ( E ) XXr, (c- 42) 



where f = 2/3Acos6>. Since the sum of & m ^ m gives U as shown in ( |C-23| ) and (C-33)-(C-34), the remaining 
problem is the sum of -B m £ m , which comes from the cluster integrals of the diagrams including two fixed 
charges. In the ring approximation, only the diagrams shown in Fig. 11 (a) are relevant. The double circle 
in the diagram denotes the sum of the diagrams as shown in Fig. 11 (b), and renormalize the fugasity £ as 

^zjTmb m r- X =^- (C-43) 
By using (C-34), it is proved that the renormalization is given by £ — ► n, where n is the particle density. 
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Fig. 11. (a) The representative diagram relevant to the optical conductivity in the ring approximation. The renormalization 
of the fugasity £ is denoted with the sum of diagram (b). The diagram (a) can be divided into the lines with I charges as 
shown in Fig. (c). 



Finally, we consider the sum of the diagrams represented by Fig. 11 (a). The cluster integrals including 
two 'fixed' charges are summed up as 
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Here, C/ denotes the sum shown in Fig. 11(c), and is obtained as 
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From ( C-42[ ) and (C-45)-( C-46| ), the equation (2.66) is derived. Similarly, A2 (r) is derive d by removing the 
last term (— l) mi in (C-44), and by replacing the electroneutrality condition (2.13) with ( 2.29| ). 

Appendix D: Low- Temperature Expansions 

In this appendix, we discuss the low-temperature expansions for the specific heat and the DC conductivity 
in the case of the ohmic damping. This expansion is obtained by asymptotic expansions in terms of the 
inverse temperature (5 = 1/T. We begin with the partition function Z expressed by 
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Here, U{6) = U{n{&)) is given by 

U(n) = n + logT(Kn + 1) - Knip(Kn + 1), 
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wh ere T (z) is the gamma function, and tp(z) is the polygamma function. At low temperatures, the integral 
in (D-l) is determined by the contribution around 9 — 0. We expand U(9) as 



17(0) ~ U(0) + ^U"{0)9 2 + 0{9 4 ), 



(D-3) 



and the integral (D-l) is replaced with the Gaussian integral by ( D-3 ), approximately. By using U"(0) = —tiq 
and asymptotic forms of T(z) and ip{z), we obtain 

Z = (1 - K)n + const. + — — + 0(njT 2 ), (D-4) 

bK Uq 

where no = n(9 = 0). We should note that the asymmetry part U""(0) also gives a term of order of 1/uq. 
However, the term takes a finite value for K — » 0, and for weak damping (K <C 1), and this contribution is 
less than the third term of (D-4). From (3.8) and (3.9), no is determined by 

2/3A = n e~ Q( " o) 



Q(n )=K(il>(Kn + l)+ f f). 
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Next, we calculate the asymptotic form of (D-5). By using A = A((3oj c /2it) K and A c g = 
A(A/lo c ) k / ( - 1 - k \ the lefthand side of (D-5) is represented as 4tt(PA cS /2tt) 1 ~ k . For (3 -> oo (n -> oo) 
the righthand side of (D-5) is expanded as 
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We substitute no = ai/3 + ao + ci-i/f3 + • • • to (D-7), and determine the constants oi, clq, cj_i so that (D-7) 
gives the lefthand side of (D-5). Thus, we obtain 
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we obtain the low-temperature expansion of the partition 



where c = {<iirKe^) K / {1 - K \ From flp"l| ) and (|D 
function, ( 3.10 ). 

In order to calculate <tdc, it is sufficient to consider the long-time behavior of S(t;8) and R(t;6) defined 
in (3.27)-(3.28). From (3.28) and JO^ ), we obtain for f -> oo 
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Here, we have used 7(6*) = 2itKn{9) / (3. From ( 2.69 ), the real-time correlation function in the long-time limit 
is obtained as 



A(t) 



Z J_ v 2tt f3 2 cos 2 ( 
Then, the DC conductivity is calculated as 
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where <7p C = e 2 a 2 /2ttK . The low-temperature expansion is obtained by expanding (D-12) over l/n(<C 1) as 
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By using the Gaussian integral formulas and an asymptotic form of the partition function 
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we obtain 
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From (D-8) and (D-15), we obtain (3.37) 
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